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Abstract 

We report the use of photoassociative spectroscopy to determine the ground state s-wave scatter- 
ing lengths for the main bosonic isotopes of strontium, ^^Sr and ^^Sr. Photoassociative transitions 
are driven with a laser red-detuned by up to 1400 GHz from the ^5*0 - ^Pi atomic resonance at 
461 nm. A minimum in the transition amplitude for '^^Sr at — 494ib5 GHz allows us to determine the 
scattering lengths 610 oq < age < 2300 ao for ^^Sr and a much smaller value of — 1 ao < ass < 13 oq 
for s^Sr. 
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Photoassociative spectroscopy (PAS) of ultracold gases, in which a laser field resonantly 
excites colliding atoms to ro- vibrational states of excited molecular potentials, is a powerful 
probe of atomic cold collisions Q]. Transition frequencies have been used to obtain dispersion 
coefficients of molecular potentials, which yield the most accurate value of the atomic excited- 



state lifetime 0, 0, 
ground _st ate atoms 
length 



^ Transition amplitudes are related to the wave-function for colliding 
I, and can be used to determine the ground state s-wave scattering 



The s-wave scattering length is a crucial parameter for determining the efficiency of 
evaporative cooling and the stability of a Bose- Einstein condensate (BEC). It also sets the 
scale for collisional frequency shifts, which can limit the accuracy and stability of atomic 
frequency standards. 

The cold collision properties of alkaline-earth atoms like strontium, calcium, and magne- 
sium, and atoms with similar electronic structure such as ytterbium, are currently the focus 
of intense study. These atoms posses narrow optical resonances that have great potential 
for optical frequency standards [3, ll^ llfil jl?! . Laser-cooling on narrow transitions is 



an efficient route to high phase-space density 



HQ, 



and a BEC was recently produced 



with ytterbium Fundamental interest in alkaline-earth atoms is also high because their 



simple molecular potentials allow accurate tests of cold collision theory 



PAS of calcium 



) jugnpecauf 



11 1 and ytterbium |lj] was recently used to determine s-wave scattering 
lengths of these atoms. This paper reports the use of PAS to determine the ground state 
s-wave scattering lengths for the main bosonic isotopes of strontium, *^Sr and ^^Sr, which 
have relative abundances of 10% and 83% respectively. We find a huge scattering length for 
^^Sr of 610 flo < (^S6 < 2300 oq. Appreciable uncertainty comes from the value of Cq for the 
ground state potential. In contrast, for ^^Sr we find — lao < ass < 13 oq. From the data, we 
also make an improved measurement ^ of the 5s5p^Pi atomic lifetime (r = 5.25 ± 0.01 ns). 
PAS results for ^®Sr, yielding a similar value of r and a different value of ass, were recently 
posted Q- 

For PAS of strontium, atoms are initially trapped in a magneto-optical trap (MOT) oper- 



l2q . We are able to produce pure 



ating on the 461 nm ^5*0 - ^Pi transition, as described in J,| 
samples of each isotope from the same atomic beam due to the intrinsic isotope selectivity 
of a MOT. For '^^Sr and ^^Sr respectively, about 7 x 10'^ and 2.5 x 10^ atoms are trapped 
and cooled to 2 mK. 
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After this stage, the 461 nm laser-coohng hght is extinguished, the field gradient is reduced 
to 0.1 G/cm, and 689 nm hght for the ^5*0 - ^Pi intercombination-hne MOT 26| is switched 
on. This MOT consists of three retro-refiected beams, each with a diameter of 2 cm and 
intensity of 400 — 800 /iW/cm^. Initially, the frequency of the laser-cooling light is detuned 
from atomic resonance by about -1.3 MHz, and spectrally broadened with a ±1.0 MHz sine- 
wave modulation. 

During a 50 ms transfer and equilibration period, the field gradient and spectral modu- 
lation are linearly ramped to 0.8 G/cm and ±0.7 MHz. The detuning is ramped to about 
-0.9 MHz. This yields about 4 x 10'' '^^Sr atoms at a temperature of 5 ± 1 fxK or 1.5 x 10^ 
^^St atoms at a temperature of about 8 ± 2 /iK. The peak density for both isotopes is about 
2 X lO^^cm"^. The intercombination-line MOT parameters are then held constant during 
an adjustable hold time. In the absence of PAS, the lifetime of atoms in the trap is about 
500 ms, limited by background gas collisions. 

The size of the atom cloud, the number of atoms, and thus the peak density are pri- 
marily determined with absorption imaging using the ^5*0 - ^Pi transition. We image along 
the direction of gravity, and transverse 1/^/e density radii are cr = 250 /im for ^^Sr and 
a = 400 /im for ^^Sr. To obtain information on the cloud dimension along gravity ^27(] . 
an additional camera monitors fiuorescence perpendicular to this direction. The cloud is 
smaller by approximately a factor of two in this axis. 

To excite photoassociative resonances, a PAS laser, tuned to the red of the atomic ^5*0 - 
^Pi transition at 461 nm, is applied to the atoms during hold times of 200—800 ms. When the 
PAS laser is tuned to a molecular resonance, photoassociation provides a loss mechanism for 
the MOT, decreasing the number of atoms. The PAS laser has negligible effect on the atom 
cloud size or temperature. When the PAS light is tuned far from a molecular resonance, it 
has no effect on the rate of loss from the trap. 

PAS light is generated from an extended cavity diode laser at 922 nm using second- 
harmonic generation in a linear enhancement cavity The laser linewidth is 80 MHz in 
the blue on a millisecond time scale, which contributes significantly to the observed PAS 
resonance linewidths. The laser frequency is measured with a wavemeter that is regularly 
calibrated with the cooling laser whose frequency is locked to the atomic ^Sq - ^Pi transition. 
The resulting 1-sigma statistical uncertainty for frequency measurements of the PAS laser 
is 200 MHz. There is also a comparable systematic uncertainty due to the sensitivity of the 
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FIG. 1: Selected region of the ®^Sr PAS spectrum. The detuning is of the PAS laser with respect 
to the atomic ^5*0 - ^Pi transition frequency. The vertical axis is the fraction of atoms remaining 
at the end of the hold time in comparison to the number in the absence of photoassociation. PAS 
laser intensity is 800 mW/cm^, and hold times range from 350 to 450ms. The minimum of the 
transition amplitude occurs when the Condon radius corresponds to a zero of the ground state 
wave function. The insets show the quality of the data and the fits to Eqs. 1 and 2. 

wavementer to the alignment of the laser into the device. 

The available PAS laser power is as high as 20 mW. Several beam geometries are used, 
but typically the PAS beam is retroreflected, with a 1/e^ intensity radius of about w = 1 
mm, yielding a maximum intensity of 3 W/cm^ on the atoms. Lower intensities are obtained 
using an AOM. For some data, a quarter-wave plate is inserted in the beam path after the 
first pass through the atoms to prevent the formation of a standing wave. This variation 
has no noticeable affect on the PAS intensity and frequency measurements. 

The observed PAS spectrum is relatively simple because the bosonic isotopes of strontium 
lack hyperfine structure. In addition, at the low temperatures of the intercombination-line 
MOT, only s-wave collisions occur so only J = 1 levels are excited. Ground state ^5*0 atoms 
collide on a potential, and, of the four states converging to the + ^Pi asymptote 
j^lj, only the state is excited in this work. Figure [T] shows spectra recorded for ^^Sr in 
a region around 500 GHz to the red of the atomic transition. The variation in transition 
amplitude will be discussed below, but first we describe the method used to quantitatively 
analyze the spectra. 

The atom density varies as n = —/?(/, f)n'^ — Fn, where / and / are the laser frequency 
and intensity. This implies the number of atoms in the trap as a function of hold time, t, 
follows 



N{t) = , (1) 



2V2r 
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where A'^o and hq are the number and peak density, respectively, at the beginning of the hold 
time. Ipk is the peak laser intensity, and ( = w"^/ (2cr^ + w"^) accounts for the laser-atom 
overlap. We approximate the density distribution as a Gaussian. The one-body loss rate 
(r) and riQ are fixed at values determined from independent measurements. Fit values of A^o 
agree well with independent measurements and are a check of the method. 

The photoassociative two-body loss rate, /3, near resonance v, is approximated as 

where /„ is the center frequency for the transition. The experimental linewidth, 7 
150 MHz, approximately equals the sum of the natural radiative linewidth for the tran- 
sition, 7„ = 61 MHz j^, and the measured laser linewidth. Including the factor 7t,/7 in 
Eq. 121 accounts for broadening beyond 71,13, 13 • ^pk is in units of mW/cm^, and is 
the collision rate constant, on resonance, for ImW/cm^ intensity from an ideal laser with 
negligible linewidth. Difficulty in accurately determining the spatial dimensions of the atom 
cloud and in aligning the PAS laser beam on the atoms leads to systematic uncertainties of 
about a factor of two in measurements of Ky (Figs. |21and|ni). Except for the most intense 
PAS transitions close to atomic resonance for ^^Sr, the linear variation with intensity used 
in Eq. |21 is a good approximation and is a constant for each transition. 

From the observed transition frequencies we can obtain an accurate value of the atomic 
5s5p ^Pi lifetime, r. For the region of the molecular potentials probed by PAS, the transition 
energies can be described with the semiclassical approximation 



E{v) = D-Xo{vd-v 
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^0 



r(4/3) 



(3) 



.2v^r(5/6)J fi^Cl 

where D is the dissociation energy, fi the reduced mass, A = 461 nm, C3 = 3^A^/167r^r, and 
f£) is a non-integer quantum number corresponding to a hypothetical level at the dissociation 
limit. Relativistic retardation, rotational energy, higher order dispersion terms, ground state 
potential curvature, thermal shifts, and all light-induced shifts are negligible at our level of 
accuracy. A fit to the ^^Sr data yields r = 5.25 ± 0.01ns. The ®^Sr data is not extensive 
enough to contribute to the determination. 

While all the data in Figs. |21 and El are valuable, the detunings of the laser corresponding to 
minima of the PAS rate are particularly important. At these detunings, the Condon radius 
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FIG. 2: Experimental and theoretical values for ^^Sr photoassociative rate constants for PAS 
laser intensities of ImW/cm^ and a narrow laser linewidth. The experimental PAS laser intensity 
used for each measurement is indicated by the symbol. Theory is for 5 fiK atoms with Cq = 
3170 a.u. (Inset) Experimental amplitudes are scaled by 1/1.4, which is reasonable given systematic 
uncertainties. The detuning of the minimum is used to determine the ground state potential as 
described in the text. The solid line is the best fit, corresponding to a minimum at -494 GHz. The 
dashed and dash-dot lines correspond to -487 GHz and -499 GHz respectively. 
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FIG. 3: Same as Fig. IS but for ^^Sr. The theoretical curve, calculated for 10 fiK atoms, is found 
using the potential determined from ^^Sr with Cq = 3170 a.u. It predicts a scattering length of 
088 = 6ao. 
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for the excitation matches a node of the ground-state wave function, which is equivalent 
to saying the overlap integral between the ground- and excited-state wavefunctions van- 
ishes. Minima positions provide precise enough information about ground-state potentials 
and wavefunctions to determine s-wave scattering lengths to high accuracy. The minima 
positions are also independent of atom density and laser intensity calibrations, and for the 
ultracold gases used here, they are independent of T, the temperature of the atoms. 

For excitation by a laser with negligible linewidth, the collision rate constant at ultra-low 
temperatures and on resonance is 5| 

""^^'''^ hQrJo (e2 + (7/2)2)' W 

where Qt = {27c fiksT/h'^)^^'^, ks is the Boltzmann constant, 7 = 7i, +7s(e, J = 0) with 
7s (e, J = 0) equal to the laser-stimulated width. At low laser intensities, 7s(e, J = 0) can 
be expressed using Fermi's golden rule as 7s(e, J = 0) = vr/rf^/eoc. Here, eo is the vacuum 
permittivity, c the vacuum speed of light and (P = \ {v\D{R)\e)\'^ where D{R) is the molecular 
dipole transition moment connecting |f) and |e), the excited vibrational wave function and 
the energy normalized ground continuum wave function respectively. Because only bound 
levels close to the potential dissociation limit are excited, 7^ is independent of v and equal 
to twice the atomic linewidth jat- By the same argument, D[R) can be approximated 
as independent of R js^ and is connected to the - ^Pi atomic dipole moment via 



D{R) = \pladat- Using dat = y^ireoh'yat^^ and the line strength factor a = y2/3 21|, we 
find D = 3.6 a.u. 

For our analysis, the inner part {R < 19 an) of ground- and excited-state potentials are 



formed by an experimental RKR potential 



33| and ab initio potential 



respectively. 



These are smoothly connected to the multipolar van der Waals expansion in C„/i?" at large 
internuclear separation. For the excited-state potential, only the C3 term contributes, and 



we use the va^ 
included IS 



ue determined above. For the ground-state potential, Cg, Cg and Cio terms are 
Relat vistic retardation effects in the asymptotic part of the excited- 



state potential are treated as described in |J]. Wave functions are calculated using a full 
quantum calculation. For the bound vibrational levels, the Mapped Fourier Grid Method 
38[ has been used whereas the ground-state continuum wave function was calculated using 
a Numerov algorithm. 

Large transition moments for ®^Sr allow us to experimentally characterize about 80 tran- 



sitions for this isotope, extending to detunings as large as -1400 GHz (Fig.|2I). This allows us 
to clearly identify minima in the transition moment at —494 ±5 GHz and —1500 GHz, corre- 
sponding to Condon radii of 62.6 ± 0.2 and 43 Oq respectively. To precisely determine the 
position of the minimum at —494 GHz, we calculated a set of potential curves with different 
node positions by varying the position of the repulsive inner-wall. Using these curves we 
performed a least-squares fit to the data shown in the inset of Fig. |21 Quoted uncertainties 
are one standard deviation. An overall amplitude factor was also varied but the best values 
were always well within the experimental uncertainties discussed above. The determined 
position of the node was independent of the value of the ground state coefficient used. 
It also did not change significantly if only data from -600 to -400 GHz was fit. 

The Sr ground state potential determined from the node position and the ground state 
Cq coefficient accurately describe the collisional properties of the system. This allows us 
to determine the s-wave scattering length age from the calculated zero-energy continuum 
wave function. However, we must carefully address the uncertainty in Ose arising from 
uncertainty in Cq. The best values in the literature are a semiempirical method yielding 
Ce = 3250 a.u. 0], and a relativistic many-body calculation of = 3170 a.u. 
although a reevaluation of 3511 inputting a more recent value of the 5s5p^Pi lifetime j2^ 
predicts Cq = 3103(7) a.u. j37|. If we take the latter value, and a minimum at —494 ± 5 
GHz, we find agg = 1450^*^370 00. (age would diverge for a minimum position of —482 GHz.) 
This large scattering length indicates a bound level in the ground state potential very near 
threshold, which implies that age is also very sensitive to Cq. For Cq = 3170 a.u. and 
Cq = 3240 a.u. we calculate age = 94:Otrmao and age = 700^90° ao respectively. The spread 
in Cq values makes it difficult to quote a rigorous statistical best value and uncertainty 
for age. Taking the extremes of the values found above, including spread in Cq and one- 
standard-deviation variation in node position, we find 610 ao < age < 2300 oq. 

PAS transitions in *^Sr are significantly weaker than in ^^Sr, and only 50 transitions were 
characterized, extending to binding energies of -174 GHz (Fig. IH)). We lack the sensitivity 
required to observe transitions to the red of the first minimum in the transition amplitude 
and thus cannot independently determine its position with high accuracy. However, the 
potential found using ^^Sr should also determine the wavefunctions for *^Sr. Because of the 
heavy mass of strontium, effects due to the breakdown of the Born-Oppenheimer approx- 
imation should be negligible as is the case for rubidium j^^. Predicted photoassociation 
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rates agree with measurements within experimental uncertainties (Fig. E)), and we calculate 
ass = 10 ± 3ao and a wave-function node at 76 ± 0.4 ao for Cq = 3103 a.u. For Cq = 3240 
a.u., Oss = 2 ± 3 ao. This yields —1 ao < ass < 13 ao. 

In conclusion, we determined the scattering lengths of '^^Sr and ^^St from PAS spectra. 
Our analysis uses full quantum calculations and is not based on semi-classical arguments 



2i 



40| that can be invalid if the scattering length becomes divergent (as in ^^Sr), small (as 
in ^^Sr), or negative. The large positive value for '^^Sr is very promising for the realization 
of a BEG through evaporative cooling in an optical trap. Such a condensate would have 
the advantage that losses due to inelastic processes would be minimal since there is no 
hyperfine interaction and spin-exhange collisions would not take place. The lifetime of a 
^^Sr condensate would be influenced mostly by three-body recombination. BEG of ^'^Sr will 
be a significant challenge |l^ because of the small scattering length, but it may allow studies 
of a very weakly interacting quantum gas with behavior similar to a hydrogen condensate 

H. 
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